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Abstract. We give a new approach to the study of statistical 



^vq mechanical systems: algebraic topology is used to investigate the 

statistical distributions of stochastic currents generated in graphs. 
O , In the adiabatic and low temperature limits we will demonstrate 

that quantization of current generation occurs. 
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1. Introduction 

In statistical physics and chemistry, especially in the study of classi- 
cal stochastic systems at the intermediate length scale, a master equa- 
tion governs the time evolution of states, in which transitions between 
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states are treated probabilistically. In its most compact form, the mas- 
ter equation is p = TpHpp, where p(£) is a one parameter family of 
probability distributions on the state space, Tz> is a constant that repre- 
sents total driving time and Hp is the master operator, which depends 
both on time t and a number j3 representing inverse temperature. 

We will be interested in varying the parameters and (3. When td 
is made large, the duration of time it takes to traverse the driving path 
is large, and one refers to this process as adiabatic (or slow) driving. 
The limiting case td — > oo is called the adiabatic driving limit. Sim- 
ilarly, one can consider the effects of low temperature on the system; 
the limiting case /3 — > oo is referred to as the low temperature limit. 

Associated with the formal solution of the master equation is an 
average current vector which represents the probability flux of a given 
initial distribution of states. In our first physics paper [CKS1], we 
argued that for generic periodic driving protocols, taking first the adi- 
abatic limit and subsequently the low temperature limit results in an 
average current vector having integer components. This quantization 
phenomenon has been observed in a variety of applications, including 
electronic turnstiles, ratchets, molecular motors and heat pumps (cf. 
the bibliography of [CKS1]). One of the purposes of the current paper 
is to give this result a mathematically rigorous foundation. Our sec- 
ond aim is to explain how algebraic topology enters the picture in an 
essential way. 

We now develop a mathematical formulation of our main results. 
Consider a particle taking a continuous time random walk on a con- 
nected finite graph T. The particle starts at a vertex i, say, and at a 
random waiting time it jumps to an adjacent vertex j where it waits 
again and so forth. Aside from the choice of inverse temperature j3, 
such a process is determined by choosing a collection of real parameters, 
one assigned to each vertex (well energies) and to each edge (barrier 
energies) of the graph. The space of these parameters is denoted by 
Air] it has the structure of a real vector space whose dimension d is 
the number of vertices plus the number of edges of T. 

Current generation occurs when the parameters are allowed to vary 
in a one parameter family. 1 We consider such a family to be parametrized 
by an interval [0,Td], in which the number represents total driving 
time. If the value of the parameters at the endpoints coincide, one ob- 
tains a periodic driving protocol; it can be represented as a pair (r D , 7) 

lr This fits with the modeling of physical and chemical processes: artificial ma- 
chines at the mesoscopic scale depend on external parameters such as electric fields, 
temperature, pressure and chemical potentials which typically vary in time. 
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in which 7: [0, 1] — > Mr is a smooth loop (equivalently, it is a smooth 
Moore loop). For each periodic driving protocol (td,j) and each (3 we 
can associate a class Q TD ,p{l) £ #i(r;R) lying m the first homology 
of the graph with real coefficients. The class is defined in terms of the 
formal solution of the master equation and is called the average cur- 
rent generated by the triple (td,j, /?). Physically, the average current 
is a measurement of the "pumping" by external forces acting on the 
system. 

The assignment 7 1— > Qrofid) describes a smooth map 

Q TD ,p: LMr -> #i(r;R) , 

where LMr is the space of smooth unbased loops in Mr with the 
Whitney C°° topology. By taking the adiabatic limit td — > 00, and 
using the Adiabatic Theorem (Corollary 13.3), we obtain a smooth 
map 

Qp: LMr -> #i(r;R) 
which does not depend on the parameter T£>. We call the latter the 
analytic current map. 

If we subsequently take the low temperature limit — > 00, it turns 
out that the resulting map is not everywhere defined. 

Definition 1.1. A loop 7 G LM. V is said to be intrinsically robust if 
there is an open neighborhood U of 7 such that the low temperature 
limit 

Q := lim Q 

is well-defined and constant on U. The subspace of LA^r consisting of 
the intrinsically robust loops is denoted by LM. T . 

The main result of this paper is a quantization result for Q. 

Theorem A (Pumping Quantization Theorem). The image of the map 

Q : LM r -> #i(r;R) 

is contained in the integral lattice Hi(T;Zj) C -£/i(r;lR). 

A version of this statement was observed earlier in our statistical 
mechanics papers [CKS1], [CKS2], and we will provide a rigorous proof 
below. A companion to the Pumping Quantization Theorem is the 
Representability Theorem, which gives a characterization of the space 
of intrinsically robust loops: 

Theorem B (Representability Theorem). There is a topological sub- 
space T> C Air such that 

LMr = L(M T \V). 
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Consequently, the space of intrinsically robust loops is a loop space. 

The subspace T> is called the discriminant, and its complement 
Air '■= M-v \ 2? is called the space of robust parameters. 

Theorem C (Discriminant Theorem). The one point compactification 
of the discriminant, i.e., T> + , has the structure of a finite regular CW 
complex of dimension dim .Mr - 2 = d — 2. In particular, the inclusion 
Air C Air is open and dense. 

Remark 1.2. A CW complex is said to be regular if its characteristic 
maps are embeddings. By [H, p. 534], such spaces have the structure of 
polyhedra. In particular, T> + is a finite polyhedron of dimension d — 2. 
We will explicitly describe the characteristic maps of T> + in §7. 

Corollary D. The inclusion LAir C LAAr is generic. In particular, 
a smooth loop 7 G LAir can always be infinitesimally perturbed to an 
intrinsically robust smooth loop 71 G LAir- 

Another main result of this paper is to give an algebraic topological 
model for the map Q: 

Theorem E (Realization Theorem). There is a weak map 

q: Air |T| 

such that the composite 

LMr -> H^Mr; Z) A #i(T; Z) 

coincides with Q, where LAir i?i(A4r;Z) is the map that sends a 
free loop to its homology class. 

(Here, |T| is the geometric realization of T. Recall that a weak map 
X — > Y is a diagram X -r- X' — > Y, in which X' — > X is a weak 
homotopy equivalence.) 

Remark 1.3. As long as T has a non-trivial cycle, the homomorphism 
g* : i?i(jVfr;Z) — > i?i(r;Z) is non-trivial (cf. Remark 7.10). In partic- 
ular, the map Q is non-trivial. 

Remark 1.4. Observe that Theorem E implies Theorem A. However, 
our actual procedure is to verify Theorem A first and thereafter use 
the tools of that proof to establish Theorem E. 

Our final result gives an interpretation of the homomorphism in 
terms of the first Chern class of a certain line bundle. Its formulation 
requires some preparation. A weak complex line bundle over a space X 
is a pair (£, h) consisting of a weak homotopy equivalence h : X' X 
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and a complex line bundle £ over X' (in terms of classifying spaces, 
this is the same thing as specifying a weak map X — > BU(1)). When 
the weak equivalence h is understood, we sometimes drop it from the 
notation and simply refer to £ as a weak complex line bundle over X. 
Since h is a cohomology isomorphism, there is no loss in considering 
the first Chern class of £ as lying in H 2 (X; Z). 

Now suppose that X = Y x Z . Then slant product with Ci(£) defines 
a homomorphism Ci(f)/: Hi(Z;Z) ->■ ff^Z) . Let S(F) = f/(l) n be 
the n-torus, where n is the first Betti number of P 

Theorem F (Chern Class Description). There exists a weak complex 
line bundle £ on the cartesian product S(T) x M. T such that 

H^Mr] Z) ^4 H\S(T); Z) = H^F; Z) 
coincides with . 
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0925618. This material is based upon work supported by the National 
Science Foundation under Grant Nos. CHE-1111350, DMS-0803363 
and DMS-1 104355. 



2. Preliminaries 
Graphs. We fix a connected finite graph 

r = (r ,r 1 ), 

where Tq is the set of vertices and Ti is the set of edges. Here we 
are allowing multiple edges between vertices and also edges linking a 
vertex to itself (loop edges). The entire structure of T is then given by 
specifying a function 

d : rwrf 

which assigns to an edge the set of vertices which it connects (T^ 
denotes the two-fold symmetric product of the set of vertices). For 
convenience, we fix a total ordering for rV Then d lifts to a map 
(d , di) : r\ — > T x T in the sense that d(e) = {d (e), di(e)}, with 
d (e) < di(e), where d (e) = di(e) if and only if e is a loop edge. The 
maps di : Ti — > T , for % — 0, 1 are called face operators. 
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The geometric realization of T is the one dimensional CW complex 
|T| given by the amalgamated union 

r u(r lX [0,1]) 

in which we identify (e, i) G T 1 x {0, 1} with rfj(e) G T for % — 0, 1. 
Populations and currents. 

Definition 2.1. The space of 'population vectors Cq(T;M.) is the real 
vector space with basis r and the space of current vectors Ci(T;R) 
is the real vector space with basis T\. If p is a population vector and 
i G r , then pj denotes the i-th component of p. Likewise, if J is a 
current vector a G T 1 then J a denotes the a-th component of J. 

The boundary operator 

d: Ci(r;R) C (T;R) 

is given on basis elements by d(a) = d (a) — di(a). Then C*(r;R) is 
the cellular chain complex of V over the vector space of real numbers. 
The spaces C;(T; M) are smooth manifolds and d is a smooth map which 
is a cellular chain analog of the divergence operator. If a current J lies 
in ifi(T;M) := ker(<9), we say that it is conserved. 

The subspace Co(T;IR) C c7o(T;IR) consisting of population vectors 
p such that J2ier P* = 1 i s called the space of normalized population 
vectors; these can be viewed as discrete probability density functions 
on the space of states r . The subspace Co(T;IR) C C7o(T;IR) of those 
p such that Pi = is called the space of zero population vectors. 

3. Driving Protocols 

Stochastic processes for periodic driving are governed by the "master 
equation" which is a certain linear first order differential equation act- 
ing on time dependent families of population vectors (see [vK, Ch. V]). 
Our master equation is a combinatorial analog of the Fokker-Planck 
equation in Langevin dynamics [vK, Chap. VIII]. 

The space of parameters. The space of parameters for T is the real 
vector space 

consisting of ordered pairs (E, W) where E: Tq — > IR and W: T\ — > W 
are real-valued functions. The function E is known as the set of well 
energies and W is is known as the set of barrier energies. We sometimes 
write (E h W a ) for the value of (E, W) at (i, a) G T x T 1 . 
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Remark 3.1. Notice that Air only depends on the number of vertices 
and edges of T, but not on the incidences. The subspace of "robust" 
parameters, which we introduce later, will depend in a crucial way on 
the incidence structure of the graph. 

Periodic Driving. A driving protocol is a smooth path 

7: [0,r D ] -^M r , 

where the real number td > plays the role of driving time. When 
7(0) = 7(td), we can view 7 as a map C TD —> Air, where C TD is the 
circle of length td. If in addition the latter map is smooth, we will say 
that 7 is periodic. When = 1, we say that 7 is normalized. 

Observe that a periodic driving protocol equivalent to specifying a 
pair 

(r D ,7) G E + x LM r 
in which 7 is a normalized periodic driving protocol. Here LAir de- 
notes the free (smooth) loop space of Air- 

The master operator. Fix a real number /3 > 0. For a given (E, W) G 
Air we can form, for each i G Tq and a G T\, the real numbers 

(1) g a = e m « , « = e m ■ 

Let g: Ci(r;M) — > Ci(r;IR) be the linear transformation given by the 
diagonal matrix whose entries are g a . Similarly, let k: Co(r;M) — > 
Cq(T; M.) be given by the diagonal matrix with entries Kj. 

Definition 3.2 (cf. [CKS1, Eq. (10)]). For a given (/?, E, W), the mas- 
ter operator is defined to be 

(2) H = -dg- 1 d*k, 

where d* : L7 (r; R) Ci(r ; E) is the formal adjoint to d. 

In particular, for fixed /3, we can view the master operator as defining 
a smooth map 

(3) H: M r ^end R (C (T;R)). 

Remark 3.3. With respect to the inner product on Co(r; M) defined by 
(u, v)k = ukv*, the master operator is self-adjoint. We infer that the 
eigenvalues of the master operator are real, and it is also easy to see 
that they are non-positive. When E = = W, the master operator is 
just the graph Laplacian —dd*. 

The master operator is also known as the Fokker-Planck operator 
to emphasize its natural interpretation as the discrete analogue of the 
Fokker-Planck operator in Langevin dynamics on smooth spaces. 
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Remark 3.4. For i,j e T , let = d~ l ({i,j}) if i ^ j and let Tj = 
{a G Til {i} C d(a)}. Setting ki jCe := g^ 1 ^, the matrix entries of the 
master operator are 



where the convention is that = when is empty, i.e., there is 
no edge connecting i and j. In particular, X^er = ^ an< ^ > ^ 
for i j (compare [vK, p. 101]). 

Remark 3.5. We offer comments on some distinctions in terminology 
between mathematics and physics. In the statistical mechanics litera- 
ture, T is usually a simple graph (no multiple edges and no loop edges). 
In this case the numbers are called rates and describe a Markov 
process on F with transition matrix H (observe that Hij = fcj jQ , with 
d(a) = {i,j} in this case). If X t denotes the state of the process at 
time t, then 

(4) H = lim P (X t+ At=j\X t = t) 

where the numerator appearing on the right denotes the conditional 
probability of transitioning to state j at time t + At, given that one is 
in state i at time t. 

Because of Eq. (1), the rates satisfy the detailed balance equation 

(5) HijKj = Hji^i . 

which states that the net flow of probability from state i to state j is 
the same as that from state j to state i. This means that the Markov 
process is time reversible [vK, p. 109]. Conversely, if the process is 
time reversible, one can show that the parameters Ki and g a are, after 
possibly rescaling, in the form given by Eq. (1). 

What we have described above is the notion of continuous time 
random walk on a graph. This is slightly more general than the notion 
of random walk considered in the mathematical literature (cf. [B, Chap. 
IX]). Mathematicians usually define a random walk to be a reversible 
Markov chain rather than the more general notion of reversible Markov 
process (the difference being that for Markov processes, one considers 
waiting times at the vertices as part of the walk). 
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The master equation. Fix a periodic driving protocol (td,j), and 
> 0. Then we have the associated one parameter family of master 
operators iJ(7(t)) £ endifj(Co(r; IR)). The master equation is given by 

(6) p(t) = r D tf ( 7 (t))p(t) • 

The master equation governs the time evolution of probability: when 
p(t) is normalized, the component p«(t) represents the probability den- 
sity of observing the state % at time t. 

The Boltzmann distribution. Suppose V is a finite dimensional real 
vector space equipped with basis B. If E: B — > M. is a function, and 
(3 > is a real number, we may form the normalized linear combination 

This is called the (normalized) Boltzmann distribution of the pair (E, f3). 
(In thermodynamics, (3 represents a multiple of inverse temperature: 
/3 = Y^fi where T is the temperature and k% is the Boltzmann con- 
stant.) The basis B identifies V with its dual space V*, so we are 
entitled to consider the function E as a vector lying in V having com- 
ponents Ei. Then for fixed (3, the Boltzmann distribution describes a 
smooth map 

p B :V^A[V], 

where A[l^] C V is the open standard simplex with respect to the basis 
T (i.e., this map sends a vector E to its Boltzmann distribution). We 
say E is non- degenerate if there is a unique j £ T such that the j-th 
component Ej of E is minimizing. 

Lemma 3.6. Let f : [0, 1] — > V be a smooth map with the property that 
f(t) is non- degenerate for every t £ [0, 1]. Then 

tends uniformly in t to the zero vector in the low temperature limit 
/3 -)■ oo. 

Proof. As [0, 1] is compact, we only need to verify the statement point- 
wise, i.e., for each t £ [0, 1]. To avoid clutter we write E{ := i% (/(£)). 
Then the z-th component of displayed derivative is 



(7) f>\ 



Case (1): i is the minimizing vertex. In this instance, the denominator 
of Eq. (7) is the square of 1 + e^ Ei ~ E ^), where each Ei — Ej < 0. 
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Hence the denominator tends to 1 in the low temperature limit. As 
for the numerator of Eq. (7), when i j, the term f3(Ej — E^)e^ Ei ~ E ^ 
tends to and when i — j it is 0. So the low temperature limit of 
Eq. (7) is 0. 

Case (2): i isn't the minimizing vertex. In this instance at least one of 
Ei — Ej is positive and Eq. (7) is dominated by k(3/e cf5 for a suitable 
choice of constants k and c with c > 0. The latter tends to zero in the 
low temperature limit by L'Hospital's rule. □ 

The Boltzmann distribution for the population space. When 
B = To, we have V = Cq(T;M.). The Boltzmann distribution in this 
case describes a smooth map 

p B : Mr -> C (T;R) 

whose value at (E, W) depends only on E and (3. It is not difficult 
to show that p B (E,W) G Co(T;IR) is in the null space of the master 
operator H(f3, E,W) (compare [vK, p. 101]). 



4. Current Generation 

For a periodic driving protocol (td,j) and (3 > 0, the instantaneous 
current at t 6 [0, 1] is defined as 

J(t) = J(/3,r D , 7 )(t) := r D g- 1 d*kp(t)eC 1 (T;R) 1 

where p(t) is the unique periodic solution of the master equation given 
by Proposition 13.1 below (here we are assuming that td is sufficiently 
large). Then the continuity equation 

dJ = -p 

is satisfied, in which J(t) plays the role of probability flux ([vK, p. 193], 
[HJ]). 

The average current generated per period is 

(8) W,r D , 7 ) := / 3{t)dt. 

Jo 

This expression measures the net flow of probability in a single period 
[0,r D ]. 
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Average current in the adiabatic limit. In the adiabatic limit 
To — > oo, both J and Q can be expressed in terms of a certain differ- 
ential Ci(T; IR)-valued 1-form A. 

For each (E, W) G Mr and (3 > 0, the negative of the restricted 
boundary map 

-d: im(g- l d*) -»■ C (T;R) 
is an isomorphism (here g~ 1 d*: C (T;R) — > Ci(T;R) and im(^ _1 9*) C 
Ci(r;R) denotes its image) Let £: C (r;R) im(^ _1 9*) denote the 
inverse transformation, and let i: im(g~ 1 d*) — > Ci(r;R) denote the 
inclusion. Then z o £ _1 defines a homomorphism 

A(P,E,W): C (r;R) ^d(r;R). 

For fixed /3 and variable (-E 1 , VF), this defines a smooth map 

A: Mr x C (r;R) -»■ Ci(r;R). 

The proof of the following is immediate. 

Lemma 4.1. T/ie map A is uniquely characterized by the following 
properties: 

(1) The composition 

M r x C {T; R) A Ci(r ; R) 4 C7 (r ; R) 

coincides with second factor projection, and 

(2) /or all J G ifi(r), we have 

(J, A), = 0, 

where (— , — zs t/ie inner product on Ci(T; R) defined by (u, v)§ = 

Remark 4.2. The operator A defines the solution to Kirchoff 's theorem 
on electrical circuits (see [B, p. 44]). Property (2) above amounts to 
Kirchoff 's voltage law with g defining the resistance matrix. 

An explicit formula for A is given as follows: choose a basepoint 
% G IV Given (W, E) G Mr and /3 > 0, define a linear transformation 
A e : Cq(T; R) — > Ci(T; R) whose value at basis elements j G r is 

(9) A e {j) = Y,Q T i J Bt J'er , 

T 

where the sum is over all spanning trees of V. The term Q^' J is the 
element of G\ (T; R) defined by the signed sum of edges along the unique 
path from % to j along T, where an edge has sign +1 if and only if its 
orientation coincides with the path. The term is the real number 
given by the T-component of the Boltzmann distribution whose vector 
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space has basis the set of spanning trees of T, where the energy function 
is given by XlaeTi W a - Then A e restricted to the subspace Co(T;M) 
coincides with A. 

Given a periodic driving protocol (td,7) and (3 > 0, application 
of the normalized Boltzmann distribution gives a loop of normalized 
population vectors 

p B ( 7 ): [0,l]->G,(r;R) 

given by 1 1 — y p B ( , ~f(t)). Taking the time derivative, we obtain a loop of 
reduced population vectors 

p B ( 7 ): [0,l]^(7 (r;R). 

Then application of A to the pair (7,p B (7)) yields a loop of currents 

A( 7 ,p B (7)): [0,l]->d(r;R) 

(This procedure describes a smooth map LM T — > LCi(r;IR).) 

The following is then a straightforward consequence of the defini- 
tions combined with the Adiabatic Theorem 13.3. 

Proposition 4.3. Let > be fixed. Then in the adiabatic limit we 
have 

lim J((3,T D , 1 )(t)=A( 1 (t),p B ( 1 )(t)) 

Tjj— >00 

and 

lim Q(/3,r D , 7 )= / A( 7 (t),p B ( 7 )(t))rft. 

By appealing to Lemma 4.1, one sees that the image of the adiabatic 
limit of Q is contained in i?i(r;R), so it defines a smooth map 

(10) Qp : LM r — ^ HxiT; R) , 

where in this notation Qp{^) ■— lim^^oo Q((3, td, 7). 

Remark 4.4. As our main results are stated in the adiabatic limit, 
there is no loss in pretending, even before taking the adiabatic limit, 
that the average current is given by the expression on the right-hand 
side of Proposition 4.3. With this change, the average current is de- 
fined without having to refer to either To or to solutions of the master 
equation. 

Definition 4.5. The map 

Q/3- LMr ->■ Hi(T; R) , 
given by Eq. (10) is called the analytical current map. 
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Fig. 1. A graph with a total ordering of its edges and its associated 
(T-spanning tree. 

5. Good parameters 

Spanning trees. For each total ordering a of the set of edges r 1; we 
may define a spanning tree T a for T by sequentially removing the links 
with the highest possible value in the ordering such that the remaining 
graph remains connected. Explicitly, let a.\ in T\ be maximal. We 
discard cti if and only if the graph Y \ a\ := (r ,ri \ {«i}) is con- 
nected. Otherwise, we retain a\. We next consider the edge a.2 which 
is maximal for r \ oc\. This is discarded if T — {a\, 0:2} is connected. 
Repeating this process, the edges which are retained form a tree T a . 

Definition 5.1. The tree T a given by the above procedure is called 
the spanning tree associated with a, or simply the a-spanning tree. 

Example 5.2. Consider the graph with total ordering a of its edges 
depicted at the top of Fig. 1. The associated cr-spanning tree is gotten 
as follows. If we remove the edge labeled 7, then the graph is connected, 
so we discard this edge. In T \ 7, the edge labeled 6 disconnects the 
graph when it is removed, so edge 6 is retained. Continuing in this 
fashion, the all edges but those labeled 3 and 7 are retained. This 
results in the spanning tree indicated by the path 1 — ?- 2 — ?- 4 — ?- 5 — ?► 6 
indicated in the bottom of Fig. 1. 

Example 5.3. A total ordering a is determined by a choice of non- 
degenerate barrier energies W: Ti — > R, where a < a' if and only if 
W a < W a i. For any such W and any spanning tree T we introduce the 
number 

w = w{T, W)= W <»- 

aeTi\T x 

When W is understood, we sometimes write w(T) for w(T, W). 

Proposition 5.4. Let W: Ti — > R be nondegenerate. Let a be the 
ordering of edges associated with W , as in Example 5. 3. Then, for any 
spanning tree T C T with T 7^ T a , there is a spanning tree T' G V , so 
thatw(T',W)>w(T,W). 
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Proof. Let (ax, . . . , and . . . , (3k) be the elements of Ti \ (T a )i 
and Ti \ T±, respectively, in the decreasing order with respect to a, i.e., 
the barrier energies are decreasing from left to right. Let j be smallest 
index such that 

• QZj 7^ f3j, and 

• cti = fa for i < j. 

Consider the spanning subgraph I"" C T, obtained from T by with- 
drawing the edges ati, . . . , <x,_i, or equivalently (3i, . . . , By the 
definition of T a , the edge aj is not a bridge of I~" (i.e., its withdrawal 
does not disconnect the graph), and Wp. < W a . for i > j. Let 
and be the two trees obtained from T by withdrawing the edge 
aj. Then there is at least one edge, say (3 S , among (3j, (3j + i, . . . , (3k that 
connects to T^ 2 \ since otherwise the edge aj would be a bridge of 
T. Therefore, by replacing the edge (3 S with aj in T results in another 
spanning tree, denoted T' that obviously satisfies the condition of the 
proposition, since W aj > Wp a . □ 

Proposition 5.4 gives an immediate characterization of a-spanning 
trees in terms of the function w (— , W). Let 7r denote the set of span- 
ning trees of V. 

Corollary 5.5. With a and W as in Proposition 5.4, the a-spanning 
tree T a is the unique maximizer of the function w(—, W) : 7r — > 

Remark 5.6. One can restate the last corollary so as to depend only 
on a: For a G Ti set W a (a) = k if a is the k-th element in the 
partial ordering given by a. Now define co(T, a) := w(T, W a ). Then by 
Corollary 5.5, T a is the unique maximizer of lo{— , a) : 7r — > 1R- 

Remark 5.7. There is a useful alternative characterizing property of the 
a-spanning tree T a associated with a nondegenerate barrier function W: 
for each withdrawn edge (i.e., an edge not in T CT ) with d(a) = {i, j}, 
the barrier W a is higher than any of the barriers associated with the 
edges of the unique embedded path which connects i to j inside T a . For 
this reason we named T a the minimal spanning tree of W in [CKS1]. 

Good parameters and the weak map q. Define an open subset 

M r C M T 

as follows: a pair (E, W) lies in Mr if and only if one of the following 
conditions hold: 

(1) there is only one absolute minimum for E : T Q — > M, or 
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(2) the function W: Ti — > R is one-to-one, i.e., the edges are dis- 
tinguished by their barrier energies. (In this instance we say 
that W is nondegenerate.) 

We call .Mr the space of good parameters. 

Let U be the set of (E, W) satisfying the first condition and let V 
be the set of (E, W) satisfying the second. Then 

M T = UU V 

where U and V are open. Each connected component of U is defined 
by specifying a vertex v E To, whereas each connected component of V 
is given by specifying a total ordering a of IY Consequently, we have 
decompositions into connected components 

U = JJ U v and V = J] K 

V cr 

For each vertex v G r , let Bi/ 3 (v) be the set of points in |T| which 
have distance < 1/3 from v in the natural metric on |T| that gives 
every edge a length of 1. Let N v C Air x |T| be the subspace given by 
U v x By 3 (v). Then the second factor projection 

N v -)• U v 

is a homotopy equivalence (it is the cartesian product of U v with 
Bi/s(v)). Now set Njj = U. V N V . Then the projection Nu — > U is 
also a homotopy equivalence. Call this projection py. 

For a given a, we let iV CT C 7W r x |T| be the subset consisting of 
V a x \T a \(l/3), where T a is the a-spanning tree and (1/3) consists of 
the points of |T| whose distance to \T a \ is < 1/3. Then the projection 
N a — > V a is a homotopy equivalence (it is the cartesian product of V a 
with a metric tree). Set Ny = U a V a . Then the projection py : Ny — > V 
is a homotopy equivalence. 

Notice that Py l {U v fl N a ) C N a . Consequently, if we set 

iV = Nu U N v , 

then a straightforward application of the gluing lemma [tD] shows that 
the first factor projection 

Px \ N ->■ 7W r 

is a homotopy equivalence. 

Definition 5.8. For good parameters, the weak map q is given by 
(11) M r ^N^\T\, 
where P2 denotes the second factor projection. 
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6. A WEAK FORM OF THE PUMPING QUANTIZATION THEOREM 

Recall the decomposition 

M r = UU V 
of the previous section, where 

u = ]lu jt v = \\v a 



3 



where j ranges through the elements of Tq and a ranges through the 
set of total orderings of V\. 

Given a loop 7 G LA4r, it will be convenient in what follows to 
think of 7 as a smooth map C — > Air, where C denotes the circle of 
radius 1/(2-71"). Let I C C be a closed arc. The contribution along / to 
the analytical current map is then given by the integral 



Jlds G Ci(iy 



where we parametrize / with respect to arc length. That is, if ii, 
is a simplicial decomposition of C into closed arcs, then 

k 



j = l Jl 3 



Assume that 7(f) C U; in this instance we say C is of type U . 

Lemma 6.1. If I is of type U, then in the low temperature limit the 
contribution along I to (^3(7) is trivial. 

Proof. On / the function E : Tq — > R has a unique absolute minimum 
v . As (3 tends to 00, the value of Boltzmann distribution p B restricted 
to arc I tends to v. This is because the component of £ G T Q in the 
Boltzmann distribution is 



and the latter tends to zero on if i 7^ v and one if £ — v when E v is the 
unique minimum. Consequently, as /3 tends to 00, the time derivative 
p B tends to zero on / (by Lemma 3.6). 

The contribution to the current of 7 along I is given by the integral 

A(7,p B ( 7 ))d S , 

(using Proposition 4.3). When (3 tends to infinity, this expression tends 
to zero. □ 
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Now consider a closed arc I = [a, b] C C with endpoints a, b such 
that 

• 7(1) C V, and 

• 7(<9J) C U. 

In this instance we say / is of type V. 

Lemma 6.2. If I is of type V, then in the low temperature limit the 
contribution along I to Qp{"f) is an element o/Ci(r;Z). 

Proof. Fix a basepoint i E T , (E, W) G Mr and (3 > 0. Recall from 
Remark 4.2 the formula 

T 

where A e : C (T;M) Ci(T;M) restricts to A on L7 (T;M). Here ^ is 
the T-component of the Boltzmann distribution for the vector space 
whose basis is the set of spanning trees of T. Recall also that the 
instantaneous current J(t) is defined as A(p B ), where p B in this case 
denotes the Boltzmann distribution for Co(r;R). Hence, inserting p B 
into the expression for A e gives 

(12) 3(t) = J2Qj'^TPl- 

T,j 

Recall that the average current is given by J(t). In particular, the 

contribution along / is given by 3(t)dt (where we are parametrizing 
/ with respect to arc length and the limits of integration come from 
the parametrization). Hence, integrating both sides of the last display, 
we obtain 

as) q=y.&' 3 r^ B - 

Since Qf^ is an integer valued 1-chain, it will suffice to prove that 

I ^ 

tends to an integer as (3 tends to oo. Using integration by parts, we 
may rewrite this expression as 

By Lemma 3.6, p B tends to zero as (3 — > oo. Hence the constribution 
to the current along / in the low temperature limit is determined by 
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the value of OrPrla- Since j(dl) C U, we deduce by the argument of 
Lemma 6.1 that the low temperature limit of OrPrla * s an hiteger. □ 

As an immediate corollary, we obtain a weak version of Theorem A: 

Theorem 6.3 (Weak Quantization). J/7 G LMr, then the low tem- 
perature limit lim^oo Qp{^) is defined and lies in the integral lattice 
Fi(T;Z) c Fi(r;R). 

7. The Discriminant Theorem and Robust Parameters 
Set 

V := Ai r \M r . 

We will first show that the one-point compactification T> + has the struc- 
ture of a regular CW complex. By definition, V is the subspace of Air 
consisting of pairs (E, W) such that E : T — > M. has more than one 
absolute minimum and W: Ti — > K is not one-to-one. 

Definition 7.1. A height function for T is a pair of functions 

ho : T -+ {1, 2} , hi:Ti^{l,...,n}, 

where 

• n > is an integer, 

• /iQ -1 (1) is non-empty, and 

• hi is surjective. 

We write h := (h , hi). 

Height functions arise in the following situation. 

Example 7.2. Given (E, W) E T>, we write ho(i) = 1 if and only if 
i is a minimum for E, and otherwise we set ho(i) = 2. We define 
hi : Ti — t- {1, . . . , n} to be the unique surjective function characterized 
by 

• hi(ai) < hi(a 2 ) if and only if W ai < W a2 , and 

• hi(ai) = hi(a 2 ) if and only if W ai = W a2 . 

The pair h := (ho, hi) is then a height function for T. 
Given a height function h = (h , hi) we define 

C(h) 

to be the set of all (E, W) G V whose associated height function is 
h, as in Example 7.2. Then V = \{ h C{h) as sets. Note that C{h) is 
non-empty if and only if h^ 1 (1) has more than one element and h^ l (k) 
has more than one element for some k. Let D(h) denote the closure of 
C(h) in V. 
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Proposition 7.3. Assume that C(h) is non-empty. Then the one- 
point compactification D{h) + is homeomorphic to a disk of dimension 
m + n, where m — 1 + |T \ 1 (1)| and n is as above. 

Proof. The space C(h) coincides with the cartesian product C (h) x 
Ci(h), where Cq{H) consists of those E : To — > R associated with ho 
and C(h\) consists of W : Tq — > R associated with h\. Likewise D{h) 
coincides with the cartesian product D(ho) x D{h\). 

Note that D(h ) is canonically homeomorphic to the space P m con- 
sisting of m-tuples or real numbers 

[X\ ) %2i ■ ■ ■ i Xrn) 

such that X\ < Xk for all A: > 1. The homeomorphism is given by 
mapping E to the map r /~ — > K, where where ~ denotes the equiv- 
alence relation on T defined by i ~ j if and only of ho(i) = 1 = ho(j). 
As r /~ has cardinality one more than the number of non-minima of 
E, such functions are identified with m-tuples of real numbers. The 
operation 

(zi, . . . , z m ) i— >■ (zi, Z2 + Zi, . . . z m + Z\) 

defines a homeomorphism to P m from the space of of m-tuples (z±, . . . , z m ) 
such that z\ e R and Zfc > for all k > 1. It is easy see that the 
one-point compactification of the latter space is homeomorphic to D m . 
Hence, D(ho) + is homeomorphic to D m . 

Similarly, the space D{h\) is canonically homeomorphic the space 
Q n consisting of n-tuples or real numbers 

{xi,x 2 , ...,x n ) 

such that x\ < X2 < ■ ■ ■ x n . The operation 

n 

(Z ± , ...,Zn)y+ (Z ± , Z! + Z2,...,y~] Zi) 

i=l 

defines a homeomorphism to Q n from the space of n-tuples (zi, . . . , z n ) 
for which z\ G R and Zi > for % > 1. The one-point compactifica- 
tion of the latter space is identified with _D n . Consequently, D(hi) + is 
homeomorphic to -D n . Finally, 

= (D(h ) x J D(/ii)) + 
= D(h ) + A D(hi) + 
= D m A D n 



□ 
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Corollary 7.4. T> + has the structure of a regular CW complex of di- 
mension d — 2, where d is the cardinality ofT II T 1 . 

Remark 7.5. With the exception of the point at oo, the open cells of 
T> + are given by the C(h) where h varies over the height functions for 
T. This is because C(h) is the interior of D(h) + . 

A top dimensional cell C(h) of T> + is given by a height function h 
in which 

• Iiq (1) has precisely two elements; 

• there is a k < n such that hi 1 (k) has precisely two elements, 
and if i ^ k we have h^ii) is a singleton for 1 < i < n. 

Robust parameters. Alexander duality applied to the inclusion T> + C 
A^jt = S d yields an isomorphism 

H d - 2 (V + ) ^ H^Mr) . 

Let C d ~ 2 (T> + ) be the cellular cochain complex of T> + over the integers 
in degree d — 2. This is the free abelian group with basis given by the 
set of (d — 2)-cells of V + . We consider the composition 

(14) C d - 2 (V + ) ->■ H d ~ 2 {V + ) = H^Mt) ^ i?i(r;R) 

where q: M. T — >• |T| is the weak map defined in Eq. (11) above. This 
homomorphism is naturally identified with a Hi(T] M)-valued chain 

0e(7 d _ 2 (£ + ;ffi(r;R)), 

and it is trivial to check that is a cycle. For a G C d ~ 2 (T> + ), let 
(a, 4>) G .£/i(r;M) denote the effect of applying the homomorphism 
(14) to a. 

Definition 7.6. A (d — 2)-cell C{h) of T> + is said to be essential if 
(C(h),(j)) G f/i(r;R) is non-trivial, where we consider C(h) as an ele- 
ment of C d ~ 2 (T> + ). A cell C(h) of T> + of any dimension is inessential 
if it is not contained in the closure of an essential (d — 2)-cell. 

Define T> to be the closure of the union of the essential (d — 2)-cells 
of V. 

Lemma 7.7. T> + is a subcomplex of T> + . 

Proof. It is enough to show that the union of any collection of top 
dimensional closed cells of T> + forms a subcomplex. Let D(h) + be any 
closed cell of V. Suppose x lies in the boundary of D{h). Then x 
lies in a unique (open) j-cell C(h'). It is straightforward to check that 
D{h') C D{h). Consequently, the boundary of D{h) + is a union of 
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Fig. 2. A two state system. 

lower dimensional cells, each of these having boundary a union of lower 
dimensional cells and so on. In particular, T> + is a union of interiors of 
certain cells, and this union is closed. Hence it is a subcomplex. □ 

Definition 7.8. Set 

M r := M r \V . 

Then we have an inclusion Air C Air- We call Xir the space of robust 
parameters. 

Example 7.9. Let T denote the graph displayed in Figure 2. In this 
case, the vector space of parameters Air is identified with IR 4 and 
the discriminant T> is identified with the diagonal inclusion M 2 C M 4 
given by (x,y) !->■ (x,x,y,y). Taking the one-point compactification of 
this inclusion, we obtain an inclusion S 2 C S A that is identified with 
T) + C Aip. The complement of this inclusion has the homotopy type 
of S 1 . Consequently, there is a homotopy equivalence 

Air ^ S 1 . 

In fact, one can make this identification precise using the loop 7 given 
by the length one periodic driving protocol ■y(t) = (cos 27it, 0, sin 2irt, 0) 
(one verifies this by showing that linking number of 7 with S 2 C S 4 is 

±1)- 

Consequently, the first homology group of Air is generated by the 
homology class [7]. The effect of g* : ifi(.Mr) _ > Hx(T) on this loop is 
to produce a homology generator of i?i(r;R) = Z (see [CKS1] for de- 
tails), so the weak map q: Air — > |T| is a weak homotopy equivalence. 
In particular, the unique 2-cell of T> + = S 2 is essential, and we infer in 
this instance Air = Air- 
Remark 7.10. With little difficulty, Example 7.9 can be generalized to 
show that q* : H\{Aiv) Hi(T) is non-trivial whenever F has non- 
trivial first Betti number. We omit the details. 

The weak map q. Let C(h) denote a (d — 2)-cell of T> + , and let 
x G C(h) be its center. We choose a small closed 2-disk meeting C{h) 
normally at x. The boundary of this disk represents a periodic driving 
protocol 7 in the space of good parameters. As an illustration, we have 
depicted such a disk in Figure 3, associated with the graph Y = XI. 
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Fig. 3. A two dimensional disk which meets a id — 2)-cell of T> trans- 
versely at its center. The green sectors are contained in U and the red 
sectors are contained in V. In this example, the current is represented 
by a non-trivial cycle so the (d — 2)-cell is essential. 



If the disk is sufficiently small then it is partitioned into four regions, 
the interior of each either contained in U or V, which is shown in Figure 
3 as green and red respectively. 2 

Each green sector corresponds to a vertex i or j of T, each giving 
a minimum for {Ek}ker - Likewise, each red sector corresponds to 
a preferred maximal spanning tree, and these are indicated next to 
each such sector together with the location of the vertices The 
topological current associated with the driving protocol 7 is defined by 
the cycle obtained by joining together two paths connecting i with j 
(cf. §8). Each path is defined by moving along one of the spanning 
trees starting at i and terminating at j (for example, in Figure 3 one 
obtains the cycle <], which is nontrivial). In particular, the cell C(h) is 
inessential if and only if this cycle is trivial. However, even more is true: 
since each path used to form the cycle is embedded, it is straightforward 
to check that the loop in T given by gluing these two paths together 
is null homotopic if and only if the two paths coincide. Consequently, 
if cr is inessential then the weak map q extends to the space given by 
attaching a two cell to .Mr along 7. If we repeat this construction 
for every inessential (d — 2)-cell, we obtain a space which is homotopy 
equivalent to the space of robust parameters Air- We have therefore 
shown that the weak map q admits an extension to the space of robust 
parameters. We denote this extension by q. 

However, in order to prove the Pumping Quantization Theorem, it 
will more convenient to give a concrete extension of the weak map q to 

2 One may justify this as follows: at x there is a unique pair of vertices i and j and 
unique pair of edges a, /3 in which Ei = Ej arc minimizing and W a = Wp. A generic 
infinitesimal perturbation into .Mr of thse values will give one of the following 
inequalities: E t > Ej, Ei < Ej, W a > Wp or W a < Wp. These inequalities 
correspond to the four sectors. 
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Fig. 4. A graph with three vertices and three edges. 

all of Air, rather than just a model for the extension up to homotopy. 
This construction will be described in the next section. 

Example 7.11. In the graph depicted by Figure 4, the edge containing 
vertices i and j is in every maximal spanning tree. Therefore, the 
discriminant T> for this graph has an inessential 4-cell. 

Definition 7.12. A barrier resolution of a height function h = (h , hi) 
is a bijection r: Ti — > {1, . . . ,\Ti\}, where \Ti\ is the cardinality of Fi, 
such that /2.1(a) < h\{(5) implies r(a) < r(/3) for all a, (3 E lY 

A barrier resolution of h enables one to associate a total ordering 
o~ r of the set I\. We saw in §5 how to obtain a spanning tree T ar for T 
associated with a r . Let 

F h = f]T ar , 

r 

where the intersection is indexed over the set of barrier resolutions r of 
h. Then Fh is a forest, i.e., a (possibly empty) disjoint union of trees. 

Proposition 7.13. A cellC(h) is inessential if and only if all elements 
°f ho (1) belong to the same connected component of Fh- 

Proof. Consider a cell C(h) that satisfies the condition that all i G 
/iq 1 (1) belong to the same connected component Th of the forest as- 
sociated with C(h). Then Th is a tree. Consider an arbitrary top 
dimensional cell C(h), so that C(h) C D(h). The cell C(h) is uniquely 
identified by two distinct vertices i, j G /i ' 1 (l), and two distinct edges 
a and with hi(a) = hi((3), together with a height function hi of 
height | Ti| — 1 with hi(a) = h^/3), i.e., W a = Wp is the only edge de- 
generacy in C(h). Then we have ^o(^) — 1 f° r k = i,j and ho(k) = 2, 
otherwise. Consider a point (E, W) G C(h), and let 7^ be a small loop 
that goes around (E, W) without intersecting C(h), i.e., staying within 
the good parameter space A^r- Using the notation of Eq. (9), we have 
for the topological current 

(15) Q(r h ) = Q> y3 '- Q? w \ 
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where h a and hp are the two possible barrier resolutions of h with 
ha, > hp and h a < hp, respectively. Then h a and hp are also barrier 
resolutions of h, and, therefore, T h C T ha and T h C , which implies 

g^C«W = gT h ,j = q^</»W Therefore, Q( 7& ) = 0, due to Eq. (15), so 
that C(h) is inessential. 

To prove the converse, let C(h) be any inessential cell. For an ar- 
bitrary pair of edges a and (3 let h% : T± — > {1, . . . , — 1} be any 
surjection such that h\(a) = hi(/3). Then there are two possible barrier 
resolutions h a and hp of h, as described above. For arbitrary distinct 
vertices i and j let ho' T — > {1,2} be the function defined by 

% ( k ) = { 1 if k = i or j, 
°^ ' ^ 2 otherwise. 

Then ft, = (ho, hi) is a height function, and C(h) is top dimensional. 
If C(h) is inessential then Qipft) = and Eq. (15) implies that the 
minimal paths that connect % to j along the spanning trees and T~ h ^ 

are identical. Furthermore, if D(h) D C(h), then C(h) is inessential. 
Consequently, the set consisting of the minimal paths that connect i 
to j inside the spanning trees associated with all barrier resolutions 
h of h consists of a single element. Denote this unique path by lh,ij- 
Then for any two distinct vertices i,j G h^^l) the path l^j belongs 
to all spanning trees T~ h , and therefore the forest associated with C(h). 
We infer that all vertices that belong to /^(l) belong to the same 
connected component of the forest. □ 



8. The weak map q 

Proposition 7.13 shows that to any inessential cell C(h) one can 
assign, in a preferred way, a tree that contains /ig 1 (1). This tree will be 
denoted Tc(h) an d referred to as the tree associated with the inessential 
cell C{h). We now cover each open cell C(h) with an open set Yh D 
C(h) that consists of all (E,W) G Air that is characterized by the 
property that if ho{i) < h (j) and hi(a) < hi(p), then Ei < Ej and 
W a < Wp, respectively, for all i,j G T and a, j5 G Ti. Setting Y = 
|J, Yh, where the union goes over all inessential cells we obtain an open 
cover 

(16) M r = UUVUY 

of the space of robust parameters. 
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For a given height function h = (ho, hi) whose cell C(h) is inessen- 
tial, we let 

N h cM T x |r| 

be the subspace given by Yh x \T C (h)\(t/3), where \T C (h)\(l/3) is the 
open regular neighborhood of \Tc\h) I given by adjoining half open subin- 
tervals of length 1/3 that correspond to those edges not in T C (h) but 
which contain a vertex of it. 

Then the projection Nh — > Yh is a homotopy equivalence. We further 
define Ny = \J h Nh, and if we set 

(17) N = Nu UN V UN Y C Mr x |r| , 
then we have a diagram 

(18) M r ^N^\T\ 

whose arrows are given by the first and second factor projections re- 
spectively. A straightforward application of the gluing lemma which 
we omit shows that the projection map pi is a homotopy equivalence. 
We infer that the Eq. (18) describes a weak map which we sometimes 
write as 

(19) q:M r ^\T\. 

By construction, the restriction of q to M. T coincides with q. 

For the sake of completeness we now sketch a proof that q coincides 
up to homotopy with the extension q of q that was described in the 
previous section. 

Lemma 8.1. The homotopy class of the weak map q coincides with the 
homotopy class of the weak map q given by gluing in 2- cells. 

Proof. The result will follow from the existence of a commutative dia- 
gram 




r 



in which the vertical maps are homotopy equivalences and the left hor- 
izontal maps (denoted p\ in each case) are also homotopy equivalences. 
The bottom maps labelled p\ and p<i comprise the weak map q. The 
top maps labelled p% and p 2 define the extension q. The space .Mr 
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is obtained by attaching suitable two cells to .Mr; it comes equipped 
with a decomposition 

M r =UUVUY 

where Y consists of the set of inessential 2-cells, each which is labelled 
as Dh, where h ranges over height functions that define a top dimen- 
sional inessential cell. Similarly, we define N C Mr x |T| to be the 
union of Nh = YhX |Tc(^)|(l/3). Then the projections onto each factor 
explicitly define the extension of q described in §7. 

The vertical maps in the diagram are given as follows. The homo- 
topy equivalence .Mr -Mr is given by the identity on .Mr and by 
mapping each closed 2-cell Dh that is attached to .Mr homeomorphi- 
cally to a small closed 2-cell Dh that intersects C(h) transversely at its 
center (the boundary of Dh is prescribed to having linking number +1 
with C(h)). A similar argument which we omit defines the homotopy 
equivalence N — > N. □ 

Remark 8.2. For the proof of the Pumping Quantization Theorem, 
it will be convenient to put the open sets Uj, V a , and Y h on equal 
notational footing. This can be done by changing the notation to 

V a :=Y h , 

where h = (h , hi) is such that ho(k) = 2 for all k G T , and h\ : Ti — > 
{1, . . . , | Til} is given by hi is defined in the obvious sense by the total 
ordering o. Similarly, we set 

U,:=Y h , 

where in this instance h = (ho, hi) is given by ho(k) = 1 for k = j 
and ho(k) = 2 otherwise, and hi is the function with constant value 1. 
Note that these notational changes necessitate a more flexible notion 
of height function, which we will call an extended height function. 

Note that the tree Tc(h) that corresponds to V a is just the a-spanning 
tree T a , whereas Tc(h) that corresponds to Uj consists of the single 
vertex j and no edges. With the above extension Eq. (17) reads N = 
N Y with N Y = U h Nh, where the union indexed over the set of extended 
height functions. 

9. The Representability Theorem 

Proof of Theorem B. It is clear that LM.r C LAir, so it suffices to 
prove the reverse inclusion. The proof will be by contradiction. Let 
7 G LM.Y- The idea is to modify 7 along a small arc in such a way 
that the value of the current changes. 
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Suppose there is an s G [0, 1] such that y = 7(5) G V. Then y is in 
the closure of an essential cell. Let e > be small. Choose a point x in 
the interior of this cell such that \y — x\ < e with respect to a choice of 
norm on M-p. Let V be an open neighborhood of 7 in LAAp on which 
Q is well-defined and constant. 

If e is sufficiently small, we can construct a smooth loop u G V 
which coincides with 7 off of (s — e, s + e), and inside this neighborhood 
oj winds once around a small disk D meeting the essential cell trans- 
versely at x in such a way that D\x C .Mr and Q(dD) is nontrivial. 
Then topologically, u is a loop obtained by concatenating 7 with dD. 
As the current Q is additive, we find that Q(u) = Q(j) +Q(dD). Con- 
sequently, Q(uj) ytz (5(7). This contradicts the assumption that Q is 
constant on V. □ 

10. The Pumping Quantization and Realization Theorems 

Consider a closed arc I — [a, b] C C of our unit length circle C such 
that 7(f) C Yfc for some h. Obviously, for c = a,b there is a well- 
defined low-temperature limit p( c > = lim^oo p B ('j(c); f3), represented 
by a normalized constant function on its support supp(p^ c - ) ) C /^(l). 
We further simplify the notation by using = Tqqi) and chose some 
arbitrary base vertex i(h) in the tree Th for each relevant height function 



Lemma 10.1. The contribution along I to the current Qp^) in the 
low-temperature limit is given by 



Proof. Using the explicit expression for the current based on the Kirch- 
hoff theorem given Eq. (12) one can show 



h. 



(20) 





TDT h 




+ 



b 



(21) 
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where Kh = {(j,T) : h (j) = 2 or Th T}. To derive Eq. (21) we first 
apply integration by parts to the explicit expression for the current, 
followed by representing the sum over the graph vertices j and spanning 
trees T as a sum over (j, T) with h (j) = 1 and T D 7\, and the 
remaining terms. We also make use of the fact that provided j G 
V'-(l), which implies j G (T h ) , and T h C T, we have Qffa = Q$$, 
i.e., the contribution to the current does not depend on the spanning 
tree T. Then Eq. (20) follows from the following properties that hold 
inside Yh, and are verified directly. For j G l {2) we have 

(22) lim pf = 0, lim dpf = 

p— >oo J B— >oo J 

and for T ^ we have 

(23) lim Qt=0, lim dg% = . 

/3— >oo B— >oo 

Since ^ T qt — 1, the properties given by Eq. (23) also imply 

(24) lim ^^ = 1, lim d V ^ = 

Eq. (20) is obtained by applying the properties given by Eqs. (22), (23), 
and (24) to the integral expression of Eq. (21). □ 

Definition 10.2. Let Q be the locally constant function given by the 
composite 

LM T -> H^Mr; Z) A if^r 1 ; Z) 
in which the first map is defined by sending a free loop to its integer 
homology class. 

Proof of Theorems A and E. Let I\, . . . , I}~ be a simplicial decomposi- 
tion of S 1 into closed arcs, with 1, . . . ,k G Z/k, and I m = [a m _i, a m ], 
so that 7(/ m ) C Y^m for some set h 1 , . . . h k of (extended) height func- 
tions. Applying Lemma 10.1, and more specifically Eq. (20), followed 
by re-grouping the terms in the sum over the arcs we obtain 



K 

(25) lim Q,( 7 ) = J2H " Q^)) 



m=l j 



(dm) 



The expression in the parenthesis on the right side of Eq. (25) does 
not depend on j. Indeed, this assertion needs only to be checked for 
for another vertex j' which lies in Th m D Th m+1 . In this instance the 
unique path running from j to j' which is contained in T^ m fl T? lm+1 

determines a one-chain c such that QJq^^ = Q^mf + c and likewise 

jryT h m+iJ' _ f^.T hm+ \,j u fyT h m,j _ ( ^ ) T hm+1 ,j _ fVr h m,j' _f- j T hm+ i,j 

Ka 'i{h m + 1 ) ~ "^ C ' n VLLm, Vj(ftm) Vj^m+l) — Vi(ftm) V^m+l) • 
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If we also account for the normalization condition for p( am \ we can 
replace summation over j by choosing any vertex j m G fl 
(h™ +1 )~ 1 (l) and then recast Eq. (25) in the form 

k 

(26) j^Q,( 7 ) = E {Q$t ~ Q^ + T) ■ 

m=l 

The right side of Eq. (26) is clearly an integer valued one-chain, so the 
proof of Theorem A is complete. 

We now turn to the proof of Theorem E. With respect to the above 
situation, consider the free loop £ : S 1 — > \T\ defined as follows: when 
the parameter seS 1 changes from the center of the arc I m to its end 
a mi £{s) goes from i(h m ) to j m along the unique minimal length path in 
the tree T h m. When s changes from a m to the center of the arc I m +i, £(s) 
goes from j m to i(h m+1 ) along the unique minimal length path in the 
tree T hm +i. It is easy to see that the right side of Eq. (26), considered as 
an element of ifi(|r|), is the image of I under the map L\T\ — > Hi(\T\) 
that associates with a free loop its corresponding homology class. On 
the other hand it is also easy to see that (7,^) : S 1 — > (Air x |T|) has 
image in C Air x |T|, and we infer (7, £) G LN. By a straightforward 
inspection of the definitions we see that the right side of Eq. (26) is 
given by Q("i) (as defined above in Definition 10.2), which completes 
the proof. □ 

11. The Chern class description 

The canonical torus. Set 

C\T;U(l)):=U(l) r > i = 0, 1 , 

where the right side denotes the set of functions Tj — > U(l). The Lie 
group 

G r := C°(r;[/(1)) 
is called the gauge group; it acts on C 1 (r; U(l)). The action is defined 
by 

(9 •/)(«) = g(do(a))f(a)g(di(a)) \ 
where g G Gr and / G Ct(T; £7(1)). Let 

^(r;£7(l)) 

denote the orbit space of this action (alternatively, let 5 : Co (T; £7(1)) — > 
d(r; £7(1)) be given by 6{g){a) = (/(do(a))</(di(a)) _1 , then H^T; £7(1)) 
is the cokernel of S). Then f£ 1 (r;£7(l)) is an n-torus where n is the 
first Betti number of T (this is the torus S(T) appearing in Theorem 
F). 
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Observe that an element of H 1 ^; U(l)) is represented by a function 
X-.Tx^Uil). 

Lemma 11.1. There is a preferred isomorphism 

H\H\Y-U{l))-Z) = H X {Y-Z). 

Proof. For a G r 1; let n a : C 1 (r; U(l)) — > U(l) denote the coordinate 
function given by restriction to {a} C I\ (use (^({a}; t7(l)) = U(l)). 
The operation a t— )■ 7c a extends linearly to an isomorphism of abelian 
groups 

d(r ; z) = [c\r- u(i)), u{i)\ = h\c\t- c/(i)) ; z) . 

We also have a similar isomorphism Co(r;Z) = H l (C°(T; £7(1)); Z). 
With these identifications, the boundary operator d: Ci(T; Z) — > Co(T; Z) 
is given by restriction 5*: H l {C l (Y; U(l)); Z) if^C ^; 17(1)); Z). 
Hence ifi(r;Z) is identified with the kernel of 5*. But the inclusion 
if 1 (if 1 (r; Z) C ker(£*) is clearly an isomorphism. □ 

A combinatorially defined line bundle. We refer the reader to 
discussion of §8, especially Remark 8.2. Recall that 

M T = {jY h 

h 

is a covering by open sets where h = (h , hi) ranges over extended 
height functions. Associated with Yh one has a tree Th : = Tc(h) such 
that Zip C Tc(h)- Fix a basepoint vertex i for Th (cf. Lemma 10.1) 

For (E, W) G Y/i and A G C 1 (r; t7(l)), we associate a complex line 
in C n , where n is the cardinality of Y\. For any vertex j of the tree Th, 
we have a minimal path p. Tw from i to j which is contained in T; this 
path defines the integer value 1-chain Qj h '^ (cf. Remark 4.2). 

Let C[r ] denote the complex vector space with basis r . Then we 
obtain a non-zero vector 

(27) v = v(h,E,W,\):= J2 (e~^ \ s ^)j G C[r ] 

je(T h ) aEPf h ' j 

where s(a) = ±1 according as to whether the direction of the path co- 
incides with the orientation of a (this sign coincides with the coefficient 
appearing of a in Q^ h '^). 

Let Cf C C[r ] denote the complex line spanned by this vector. 

Lemma 11.2. If we choose a different basepoint vertex the complex 
line Cf remains unchanged. 
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Proof. The amalgamation of the minimal length paths pj h ' 1 and Pj h '^ 
produces a new path pT h,i pj w from i to j. If this path is minimal, 
then it is Pf h ^ and clearly, we have 

n a « q) =( n A « a) )( n xs ^- 

a &pJ h ' 3 aeP^' 1 ' «ePj hJ 

If the amalgamated path isn't minimal, then this formula still holds 
because the factors corresponding to indices occurring more than once 
cancel. This gives independence with respect to the basepoint vertex, 
as the first factor on the right is independent of j. □ 

For fixed h the assignment (E, W, A) i-> Cv(h, E, W, A) describes a 
line bundle ^ over C l (T; U(l)) x Y^. In fact, it is straightforward to 
check that ^ is trivializable. We now use the clutching construction 
to glue these line bundles together as h varies. This will produce a 
line bundle over £ over C^fT; U(l)) x .Mr- To check this, it suffices to 
establish the following. 

Lemma 11.3. Given height functions h and h! , let 

a: C\V- £/(l)) x (Y h n Y w ) -> C\T; U(l)) x Y h 

and 

b: C\T; U(l)) x (Y h fl K h -> C 1 ^; 17(1)) x 
denote the inclusions. Then there is an isomorphism of line bundles 
4>ab'- b*^h' a*£h- Furthermore, this isomorphism satisfies the cocycle 
condition (p ac = 4> a b4>bc- 

Proof. Associated with (X,E,W) G C 1 (r;f/(1)) x Yh and a basepoint 
vertex i for Th H TV, we have a non-zero vector v which is defined 
by Eq. (27). To indicate the dependence of this vector on h, let us 
redenote it by Vh- Similarly, for (A, E, W) G C l {T; U (1)) x Yh we have 
Vh>. Then define 4> a b(z 4 v h') '■= z • v^. The cocycle condition is then 
immediate. □ 

Let the gauge group Gr act diagonally C 1 (r; Z7 (1)) x Air (where the 
action on the second factor is trivial). Then Gr also acts in an evident 
way on the total space E(£) of the line bundle £ equipping it with the 
structure of a Gr-equivariant line bundle. Taking orbit spaces defines a 
line bundle £ over 7f 1 (r;t/(l))x Air- Ktt: C 1 {Y- 1 U{1)) -> /f^r^l)) 
is the quotient map, then £ is given by the base change of £ along 

7T x id: C^T; 17(1)) x Al r -> H^T; (U(l)) x Af r • 
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Naturality of Chern classes gives a commutative diagram 



(28) 



#i(.M r ;Z) 



ci(0/ 



>Ci(r ; z) 




TT 



Here we have used the preferred isomorphism if 1 (iJ 1 (r; U(l)); Z) = 
Hi(T;Z) of Lemma 11.1 as well as a similarly constructed identifica- 
tion J ff 1 (C 1 (r;f/(l));Z) = Ci(T;Z). With respect to these identifica- 
tions, tt*, which is the homomorphism induced by it on first integer 
cohomology, is just the canonical inclusion Hi(T; Z) C Ci(T; Z). 

Theorem 11.4. The homomorphism ci(£)/ coincides with q*. 

In order to prove Theorem 11.4 we digress to explain how holonomy 
relates to the homomorphism given by slant product with the first 
Chern class. If £ is a complex line bundle over a connected space B 
and structure group U(l), we have the holonomy map 



In fact, hj: can be chosen in such a way that if we choose a basepoint of 
B and restrict to the based loop space QB, we can deloop to a map 
B — > BU(1) = CP°° that classifies the bundle £ and hence the Chern 
class ci(£). Consequently, if we choose fi£ in this way, it determines the 
Chern class. 

Now suppose B = X x Y . Then we can restrict to the subspace 
X x LY C LX x LY = L(X x Y) and take the adjoint to obtain a 
map 



where F(X, U(l)) is the function space of maps from X to U(l). Then 
the diagram 



(29) 



hz'.LB-* U{1) . 



LY — >■ F(X,U(1)) , 



(30) 



LY 



>F(X,U(1)) 



H X {Y-Z) 



> H 1 



(X;Z) 



commutes, where the left vertical map sends a loop to its homology 
class and the right vertical map sends a function to its homotopy class 
considered as an element of i? 1 (X; Z) = [X, U (1)]. 
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Proof of Theorem 11. 4- By Eq. (28) it suffices to show that Ci(£)/ co- 
incides with the homomorphism 

H x (Mv] Z) A Hx(T; Z) ^ d(T; Z) . 

Suppose we are given (A, 7) G C 1 (r; U (1)) x LM.Y- As done previously, 
we partition S 1 into closed arcs [a^, dk+i] for < k < n with a n+ i = a , 
such that the projection of such an arc into .Mr is contained in a 
neighborhood of type Yj lk _ 1 . Choose a basepoint vertex lying in the 
intersection T& RTfc+i, where Tk denotes the tree associated with Yh k . 
Then we have a minimal length path P^ fe,ifc+1 from to ik+i and the 
product 

n 

(31) n n K (a) eu(i) 

k=0 ^p^.ifc+l 

describes a map C l (T;U(l)) — > f/(l) that gives the holonomy around 
7 (here we are using Lemmas 11.2 and 11.3). 

In the special instance of A G C 1 (r;[/(1)) which is identically 1 
except for a single edge a, the value of the map C 1 (r; U(l)) — > U(l) 
at A is given by A* where q represents the net number of times a is 
traversed, with orientation taken into account. Consequently, if we 
identify H^C 1 ^; U(l)); Z) with d(T;Z), then A« is identified with 
the chain qa. It follows that the map C 1 (r; U(l)) — > U(l) defined by 
Eq. (31) corresponds to the integer cycle in Ci(T; Z) given by 

n n 

(32) £ £ S (a)a := 

fc=0 agp T fc ,i fc+1 fc=o 

On the other hand, the paths Pj k k ' lk+1 describe a lift of 7: S l — > M r 
through the space N appearing Eq. (17) (roughly, one defines the lift by 
mapping the midpoint a' k of the arc [a*., afc+i] to the point {y(a'k),ik) 
and uses p^ k ' tk+1 to connect these points). Then application of the 
projection map P2 : N — > \T\ to the given lift produces a map S 1 — > \T\ 
that represents 9* ([7]) G Hi(T;Z) C Ci(T;Z). From this description 
it is straightforward to check that 9* ([7]) coincides with the element 
defined by Eq. (32) (cf. Eq. (9)). □ 

12. The ground state bundle: a conjecture 

By coupling the master operator with elements of the torus H 1 (T; U(l)) 
one can extend the master operator to a self-adjoint operator over the 
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complex numbers. This extension is called the twisted master opera- 
tor; its eigenvalues are real and non-positive. The eigenspace associated 
with the maximum non-zero eigenvalue is called the ground state. One 
may use the twisted master operator to define another weak complex 
line bundle 77, this time over iJ 1 (r;f/(l)) x .Mr, where Air C Air 
is characterized by the condition that the ground state at each point 
is non-degenerate, meaning that it has rank one. Then roughly, 77 is 
defined by taking the ground state at each point of the base. We call 
this the ground state bundle. Arguments from physics suggest that 
the ground state bundle is equivalent to the weak complex line bundle 
£ that was defined in the previous section. In what follows we will 
formulate this idea as a pair of conjectures. 

The twisted master operator. The twisted master operator, defined 
below, is a smooth map 

H: C^T; 17(1)) x Mr ^ end c (C (X; C)) , 

where C (X; C) is the complex vector space with basis T . It extends 
the master operator in the sense that 

H(1,0,E,W)=H(0,E,W) 

where 1 G C 1 (r; U(l)) is the function with constant value 1 G U(l) and 
where we are interpreting the right side of this identity using extension 
of scalars. 

For A G C\r; 17(1)), let A: d(T; C) ->■ d(r ; C) be given by rescal- 
ing each basis element a by X(a)a. Then 

H(X, f3, E, W) := -dg-^Xd*^. 

It is clear from the definition that H is self-adjoint. In particular, its 
eigenvalues are all real. 

Let the gauge group Gr act on endc(Co(X; C)) via conjugation and 
trivially on Air- The following is then a formal consequence of the 
definitions. 

Lemma 12.1 (Gauge Symmetry). The twisted master operator is Gr- 
equivariant, i.e., for h G Gr, we have 

H(h ■ X,/3,E,W) = h- H(X, f3, E, W) . 

In particular, for each (E, W) G Air, the spectrum of H(/3, E, W) = 
H(l, (3, E, W) is invariant with respect to the action of the gauge group. 

Definition 12.2. Let A: V — > V be a self-adjoint linear transforma- 
tion of a finite dimensional complex vector space V, all of whose eigen- 
values are non-positive. Then the ground state L is the eigenspace for 



QUANTIZATION OF FLUCTUATING CURRENTS 35 

maximal eigenvalue of A. We say that A has nondegenerate its ground 
state L has rank one. 

An analytically defined weak line bundle. Define an open subset 

Mr C IR+ x Mr 

to be the set of those ((3 Q ,E, W) such that for every A G C l (T;U(l)) 
and every > /3q, the twisted master operator 

H(X,f3,E,W): C (T;C) C (r;C) 

has a non-degenerate ground state. 

For (A, fto, E, W) G C 1 {Y;U{1)) x Mr, let us denote the ground 
state of the twisted master operator by L(/3 , A, E, W); it is a complex 
line in Co(T; C). Consider 

E — {(A, /3 , W, G L(/3 , A, £, W)} 

which is topologized as a subspace of C 1 (r;?7(l)) x Mr x Co(T;C). 
Then we have an evident projection 

p: E -»■ C^rjC/^)) x 7W r . 

Lemma 12.3. TTie map p is a smooth complex line bundle projection. 

Proof. Let V = Co(T; C) and let L l {V, V) denote the space consisting 
of complex linear self-maps of V having corank one. Then L l (V, V) is 
a smooth manifold of real dimension 2|r | 2 — 2 (see [GG, prop. 5.3]). 
The operation which sends a complex linear self-map to its null space 
defines a smooth map 

l\v,v) ^P x (y) 

whose target is the projective space of complex lines in V. The com- 
position 

C 1 ^- C) x Mr A L\V, V) ^ F\V) 

is therefore smooth and the pullback of the tautological line bundle 
over P 1 (V r ) gives the projection p. □ 

Let 7r : Mr -> Mr be given by the projection (/3 , E, W) (E, W). 

Conjecture 12.4. The image of it is Mr, and it: Mr —> Mr is a 

weak homotopy equivalence. 

Let f\ be the complex line bundle defined by Lemma 12.3. The gauge 
group Gr acts on both the total and base spaces making fj into a Gr- 
equivariant complex line bundle. Taking orbits, we obtain a complex 
line bundle -q over H l (T\ U{1)) x Mr- 
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Let h: H\T; U{1)) x Mr -> H^T; 17(1)) x Mr be given by id x tt. 
Then using Lemma 12.4, the pair (77, h) is a weak complex line bundle 
over H x (r;U(l)) x M r - 

Then slant product with the first Chern class of rj gives a homomor- 
phism 

cifo)/: ^(^^^^(rjZ), 

where we have implicitly used the identification ^(.MrS ^) — ^i(«A^rj ^) 
of Lemma 12.4 and also the identification H 1 (H l (T;U{l));Z) = #i(T;Z) 
of Lemma 11.1. 

Conjecture 12.5. The homomorphisms c\{rj) / andc\(£)/ coincide. 

13. Appendix: an adiabatic theorem 

Here we formulate and prove an adiabatic theorem for periodic driv- 
ing. Roughly, it states that for slow enough driving a periodic solution 
of the master equation exists and is unique, and furthermore, in the 
adiabatic limit this solution will tend to the Boltzmann distribution 
taken at the associated normalized driving protocol. 

Let us introduce the evolution operator XJ it, to) = U(t,t ; H,td) for 
< to < t < 1, which is the unique solution to the initial value problem 

-U(t, t ) = r D U(t, to)#(7(*)), U(t , t ) = I, 
at 

where I denotes the identity operator. We remark that U (t, t ) is also 
called the path-ordered exponential and is sometimes expressed in the 
notation 

f exp(r D f dt'H^t'))) 

J to 

(cf. [L]). 

Then it is elementary to show that the master equation 

p(t) = r D H(j{t))p{t) 

has formal solution 

p(t) = tf(t,0)p(0). 

Proposition 13.1. Let (tjj,j) be a periodic driving protocol. Then 
there is positive real number r such that if To > r , then there is a 
unique periodic solution p(t) to the master equation, i.e., p(0) = p(l). 

Proof. We shall use abbreviated notation and write p B {t) in place of 
p B (7(£)). For any solution to the master equation pit), set 

m := P B (t) - P(t) . 
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Then £: [0,1] —> (%(r;R) is a family of reduced population vectors. 
Furthermore, is periodic if and only p(£) is, and 

(33) P(t) = P B (t)+m- 

Inserting Eq. (33) expression into the master equation and using the 
fact that the Boltzmann distribution lies in the null space of the master 
operator, we obtain the first order linear differential equation in £, 

(34) i(t)-T D H(t)£(t) = -fP(t), 

where H(t) is shorthand for H(j(t)). 
Applying U(l,t) to Eq. (34) we get 

U{l,t)£-U{l,t)T D HZ=-U(l,t)fP. 

Notice that the left side of the last display is just d/dt(U(l,t)£). Inte- 
grating both sides we obtain 

U(l,t)£(t) = - f U(l,t')p B dt' + C . 
Jo 

Setting t = we see that U(l, 0)£(0) = C. Evaluating at t = 1 and 
using the fact that £7(1, 1) = I yields 

£(l)-tf(l,0)£(0) = - [ U(l,t')p B dt'. 

Jo 

Consequently, ^(0) = £(1) if and only if 

(35) (/ - 17(1, 0))f (0) = - / £7(1, t')p B dt' . 

Jo 

It is therefore sufficient to show that the operator I — U(l, 0) is invert- 
ible, when considered as an operator acting on the invariant subspace 
Co(r;R), provided td is sufficiently large. 

Let A and c be the constants obtained in Lemma 13.2 below. If we 
set t := (1/A)ln(2c), then we have ||£/(1,0)|| < 1/2. It follows that 
I - U(l, 0) is invertible on C (T; R). □ 

The last part of the proof of Proposition 13.1 rested on an estimate 
that appears below. To formulate it we use the norm on Co(r; R) given 
by || £|| = ■sj (£, £) where the inner product is the one induced by the 
standard inner product on Co(r;R). If A is an operator on Co(r;R) 
then we define ||A|| := sup^ ||||f H" 1 = sup €=1 ||A£||. 
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Lemma 13.2. For a periodic driving protocol (td, 7), there are positive 
constants X and c such that for all t, t e [0, 1] we have 

\\U(t,t )\\ <ce -^(*-*o). 

Proof. Consider the time-dependent inner product k% = «(7(t)) in 
C (T;R), denned by K t (£, V ) = E i6 r e pE ^%Vj (the is just (Z, V ) kt 
in the notation of Remark 3.3). Then for all t the operator H{t) = 
if (7(i)), when considered as acting on (%(r;R), is self-adjoint with 
respect to the inner product K t and its spectrum is strictly negative. 

Set A := — sup t6 r 01 ] <r(H(t)), where o~{T) denotes the spectrum of 
a linear operator T. Then A > 0, and the spectrum of the opera- 
tor H (t) = H{t) + XI is non-negative for all t. Let £/o(t, to) be the 
corresponding evolution operator. Then U(t,t ) = e~ ATD< * _ ' )£7o(t, t ). 
Hence, \\U(t,t )\\ = e- ATo( *- io) ||£/ (Mo)||- So all we need to prove is 
that ||{7o(£, to) || is uniformly bounded. 

Let r](t) be the solution of the master equation rj(t) = T D H (t)r](t) 
with the initial condition 77 (to) = £• We then have 

j^rfit),^)) = k t (r ] (t),r ] (t)) + 2r D K t (H (t) V (t),r ] (t)) 

(36) < kt(ri(t),Ti(t)), 

since for all t we have K t (H (t)r)(t),r)(t)) < 0. 

Since 77 (t) 7^ provided r)(t ) = £ 7^ 0, we infer that K t (r](t),r](t)) > 
for all t. By compactness, \\n t \\ is bounded below, and since \\kt\\ 
is bounded above, we infer that there is a constant A > 0, so that 
Kt(v(t), 7 /(t))(^t(r/(t), r/(t))) _1 < A. Combined with Eq. (36) this im- 
plies (d/dt) In K t (rj(t),r](t)) < A, and further implies the uniform bound 

/ 37 x K t (U (t,t )£,U (t,t )£) = Kt(v(t),y(t)) < e A(t-t ) 

The uniform bound provided by Eq. (37) implies the uniform bound 

(38) M) < B , 

with respect to the standard inner product for some B > 0, which 
immediately implies the uniform bound ||C/o(t, to)|| < B. □ □ 

Corollary 13.3 (Adiabatic Theorem, cf. [vK, V.3]). Let (td,7) be a 
periodic driving protocol, with sufficiently large. If p(t) denotes the 
periodic solution of the master equation, then 

p B (l(t))= lim p(t). 

T£>— >0O 
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Proof. It is enough to show that lim Tl5 _ s . 00 £(£) = where is as in 
the proof of Proposition 13.1. We first show that lim^^oo £(0) = 0. 
To see this, start with the estimate 



U(t,t')p B (t')dt'\\< / \\U(t,t')\\\\p B (t')\\dt' 



<c sup ||p B (t')ll / 
t'efo.il Jo 



-*T D (l-t) dt 



< 



'£[0,1] 
ac 



where a = sup t / 6 r 01 i ||p B (t')ll- Recalling that ||£7(1,0)|| < 1/2, we have 
|| (7 - £7(1, 0))- 1 1| < 2. Consequently, 

||£(0)|| = ||(J-17(1,0))- 1 fu(l, t ')fPdt'\\ 

Jo 

< \\{i-u{i^)y l \\\\ f u(i,t')p B dt'\\ 

Jo 

lac 

< 



Therefore, Hm T£) _ ) . 0O £(0) = 0. 

The proof that lim^^oo £(t Q ) = for any t G [0,1] is similar, 
using a suitable modification of the above estimate with t i n place of 
and Lemma 13.2) to give a similar bound for ||£(£o)|| (we omit the 
details). □ 
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